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Abstract 

The governed equations for the order parameter, one-time and two-time correlators are ob- 
tained on the basis of the Langevin equation with the white multiplicative noise which amplitude 
is determined by an exponent < a < 1 (x being a stochastic variable). It turns out that 
equation for autocorrelator includes an anomalous average of the power-law function with the 
fractional exponent 2a. Determination of this average for the stochastic system with a self- 
similar phase space is performed. It is shown that at a > 1/2, when the system is disordered, 
the correlator behaves non-monotonically in the course of time, whereas the autocorrelator is 
increased monotonically. At a < 1/2 the phase portrait of the system evolution divides into two 
domains: at small initial values of the order parameter, the system evolves to a disordered state, 
as above; within the ordered domain it is attracted to the point having the finite values of the 
autocorrelator and order parameter. The long-time asymptotes are defined to show that, within 
the disordered domain, the autocorrelator decays hyperbolically and the order parameter be- 
haves as the power-law function with fractional exponent —2(1 — a). Correspondingly, within the 
ordered domain, the behavior of both dependencies is exponential with an index proportional 
to —tint. 
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1 Introduction 



Within the framework of the ordinary thermodynamic approach, it is postulated that the 
bath is passive with respect to a variation of the system state parameter x which is an 
amphtude of the hydrodynamic mode [1]. In such a case the noise of the corresponding 
stochastic process x{t) is additive one in the course of time t to have the temperature as 
intensity, which is independent on the variable x. On the contrary, the amplitude g{x) of 
a multiplicative noise varies with the stochastic variable x. Such type examples represent 
a population dynamics [2], directed percolation [3], Levy flights [4] and so on. 

According to Ref. [5], the multiplicative function g{x,t) represents the homogeneous func- 
tion for the systems with a self-similar phase space. In such a case we can write the noise 
amplitude as the power-law function 



that is singular in character because g = at x = 0. This kind assumption allows us to 
consider such models as ordinary thermodynamic system with an additive noise (a = 0), 
directed percolation process (a = 1/2), population dynamics (a = 1) and so on. 

This work is organized as follows. In Sec. 2 the governed equations for the first moment 
(order parameter) of the stochastic variable x{t) as well as the one-time and two-time 
correlators are obtained on the basis of the Langevin equation with the white multiplica- 
tive noise. It turns out that equation for autocorrelator gains an anomalous average of 
the squared power function (Q) with fractional exponent 2a (the fractional average). Sec. 3 
deals with the determination of this average for the stochastic system with a self-similar 
phase plane, whose distribution function is a homogeneous one and can be approximated 
by a power function [5]. Within the framework of this approach, the system behavior 
is governed by the value of the exponent a in Eq.(0) [5,6]. At a > 1/2 the system is 
disordered to be represented by correlator and autocorrelator (Sec. 4). At a < 1/2 the evo- 
lution of the order parameter ought to take into account (Sec. 5). Sec. 6 contains a short 
conclusions. 



2 Basic equations 

As usually, let us start with the Langevin equation 



for a stochastic variable x{t). In right-hand side of Eq.(^, f{x) is the deterministic evolu- 
tion force, the second term defines the multiplicative noise with the amplitude g{x). The 




a G [0, 1] 





(2) 
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statistical properties of the Langevin force ^ (t) are standard: 



m) = o, mat')) = 5it-t') (3) 

where the angle brackets denote the averaging. In order to study the simplest system prop- 
erties, we consider the order parameter rj(t) = {x(t)), the two-time correlator G(t,t') = 
{6x{t)Sx(t')) , Sx(t) = x{t) — {x(t)) and the autocorrelator (structure factor) S(t) = 
((5x(t))^). Within the framework of the white-noise approximation that is expressed by 
Eqs.(l^) and the Ito calculus, we can treat g{x) and ^(t) as statistically independent func- 
tions. Then, at averaging of Eq.(|D we can set {g{x)^{t)) = so that the evolution of 
the first moment is defined by the equation rj = {f{x)) where {f{x)) ^ f{v)j dot stands 
for the derivative with respect to the time t. Without loss of generality, the deterministic 
part of the evolution force can be chosen in the polynomial Landau form: 

/(^) = -^, Vix) = -'-x^ + \x^ (4) 

where —e is an external driven parameter type of the dimensionless temperature counted 
from a critical value. Performing the averaging under the force definition (^), we get term 
{x^) that is reduced to 77(0;^) = rj^r]'^ + S), in accordance with the cumulant expansion. 
As a result, the evolution equation for order parameter takes the form 

f] = ri{e — rf^ — SrjS (5) 

where the values rj, S are dependent on the time t. 

By analogy, accounting the equation 

{gixit))ximt)) = {9ixit))xit'))m) + {gixmmit') + {gixm{4mt)) = o 

that follows from Eq.(^) and above-mentioned property of the white noise, we obtain the 
equation for two-time correlator 

t') = {e- 3[v\t) + SmG{t, t'). (6) 

The problem lies now in obtaining an evolution equation for the variance S = (x^) — 77^. 
It can be performed if we use the relation dx^ = (x + dx)^ — x^ where, in accordance with 
Eq.(^, the differential dx is written as follows: 

dx = /(x)dt + g{x)dw, dw = ^(t)dt, {dwf = dt. (7) 

Then, up to the first derivative, the equation for (x^) takes the form 

A(a;2) = 2(x/(x)) + (/(x)). (8) 
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Here the last term is the average intensity of the multiphcative noise being a result of the 
interaction of the variable x with a bath, whose variables have been appropriately elim- 
inated. Inserting definition given by Eq.(0) into Eq.(Q), we come up against the average 
(x^") of the stochastic variable with a fraction exponent 2a. Further, we are coming to the 
obtaining an expression for such fractional average in terms of the cumulant expansion. 



3 Calculation of the fractional average 



Let us admit that, apart from the initial distribution P{x), there is another distribution 
Pq{x) specified by a positive parameter g < 1. Moreover, for P{x) and Pq{y) we assume 
the following relation: 

x'^P{x)dx = yPq{y)dy. (9) 

Then, the probability density Pq{y) for the new stochastic variable y = x'^ can be rewritten 
as the normalized distribution 

PM = q-V'-'^^'Piy'^')- (10) 

Denoting the averaging over the distribution Pq{y) as (• • .)g and using the designation 
(. . .) for the initial distribution P{x), we obtain 



{x') = {y)q; (11) 

(x") = / x'^P{x)dx, {y)q = q-' f yy^'-''^/'^P{y'/'')dy. 



Thus, the distribution (ITUp allows us to use the usual cumulant expansion for any average 



(x*^) with the fractional exponent q. 

For self-similar stochastic systems the distribution function can be written in a power-law 
form 

P(x)~Aa;-2^ A = ^|l-2a| (12) 

to be a homogeneous function [5], [7]. Here factor A is responsible for a cut-off procedure 
in normalization condition 

1/6 

2 J P{x)dx = 1 (13) 

b 

with cut-off parameter 6 — > 0. The integrating with the distribution function ( [T^ ) gives 
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)=aJ x"«x-2"dx = A{1 - 2a + nqy^x^-^''+'"^, (14) 
(x") = A / x"x-^"dx = A(l - 2a + (15) 



As a result, using Eq.(|TT|), we obtain the relation 

= a„(g)(x")P"(«) (16) 
where the exponent Pn{(l) and the multiplier Q;„(g) are introduced as follows: 

Pn{q) = ^ + a^q) = p~\q) (1 - 2a + n)^"(^)-^ (17) 

1 — 2a + n 

Now, we are ready to formulate the equation for the autocorrelator S = (x^) — rj'^ on 
the basis of Eqs.(|), ([l|), (p!6D, (p!?!). According to [6] a keypoint of the system with the 



multiplicative noise is that its behavior is governed by the magnitude of the exponent a in 
Eq.(|I]). Atl/2<a<l, when the fractal dimension of the phase space D = 2(1 — a) is less 
than 1, the system is always disordered and its evolution is represented by the correlator 
G{t,t') and the structure factor S(t). The former is governed by Eq.(||), for the latter it 
follows 



S = 2Si£-3S) + a2SP\ (18) 
«2 = «2(a) = p2 = P2{a) = (3 - 2ay' 

from Eqs.(|D, (§), ([l|), (0), ([T7| ) if we put q = a, n = 2. Within another domain < a < 
1/2, the above fractal dimension D > 1 so that the system can be ordered and instead of 
Eq.(p!8D we obtain 



S = 2S 



(19) 



ai ^ ai(2a) = yl^i-^a)^,^-^,^ ^ ^^^^^ ^ ^2(1 - a)]-' 

at g = 2a, n = 1 in Eqs.(|lD, (^). 



4 Evolution of disordered system 



As pointed out above, in the case when the exponent a > 1/2 (the fractal dimension 
D < 1), the system is governed by Eqs.(|l8D, @ for the one-time and two-time correlators 
S{t), G{t,t') being the structure factor and Green response function. The form of the 
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time-dependence for the former is shown in Fig. la. It is seen that S{t) monotonically 
increases to the stationary magnitude 5*0 determined by the equation 

£ - 3^0 + («2/2)^r"' = 0. (20) 
In the limit 5 < 1 when > 5 > S^, Eq.(|l8D gives the power-law time dependence 

/ A2{l~a) \P2/(1-P2) 

where we put S(t = 0) = 0. In opposite case So — S <^ Sq one has the exponential 
dependence S — Sq oc e~^^, A = 6(2 — p2)So — 2(1 — p2)£- According to Eq.(pOD the 
stationary value 5*0 raises with e increase from the minimal magnitude (0:2/6) (see 
Fig.lb). 

The solutions of Eq.(^ for different values a and e are shown in Fig.2. We plot corre- 
spondent dependencies of the Green function G{t,0) at identical initial conditions. It is 
interesting to observe that the correlator G{t,0) reaches firstly its maximum and then 
monotonically decreases to zero. The function G{t,0) attains its maximum more sharply 
if the parameter e increases (cf. curves 1, 2). The exponent a increasing drives to the 
similar effect (cf. curves 1, 3). 



5 Evolution of ordering system 



Now one has a < 1/2, D > 1 and the system behavior is governed by Eqs.(^, (^, 
([TP|), from which the first and third state the enclosed system of differential equation. To 
analyze the latter, it is convenient to use the phase plane method. As is seen from the 
corresponding phase portrait in Fig. 3a, at small values e there is only one attractive point 
?7o = 0, 5*0 = e/3 (Fig. 3a). With e increase at the point 



^0 = 7: 

2-pi 



-(2 -pi)ai 



2/(4-pi) 



f22) 



a bifurcation creates new saddle and attractive points (see Fig.4) with coordinates rjc = 
[(2 - pi)(4 - ^»l)"^eo]^^^ Sc = (2/3)(4 - pi)~^£:o- It is seen from Fig.5 that bifurcation 
temperature Eq is increased infinitely with the exponent a growth to the critical value 
a = 1/2. The coordinates of the new stationary points are determined by equations 

e-r^l- (3/4)air/g-2 = 0, So = {^a,r\e - 3SoY^'^/'^'\ (23) 



that is obtained from Eqs.(^, (|T9D at ?) = 0, 5 = 0. The corresponding dependencies 
rjo{e), So{e) are depicted in Fig.4 where the dashed curves are respective for the saddle S 
and solid ones - for the attractive point C. 
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It is interesting to note that the system undergoes the phase transition of the first order, 
despite of the bare x^-potential (Q) corresponds to the continuous one. Thus, at small 
values of the exponent of the multiplicative noise (|1|) (a < 1/2) the fluctuations transform 
order of the phase transition, whereas at a > 1/2 ones suppress the ordering process at 
all. 

Let us return now to analysis of the time dependencies of the main averages under con- 
sideration. Firstly, we analyze the system evolution to the disordered state (a vicinity of 
the point Co in Fig. 3) for extremely large time t ^ oo. The keypoint of our consideration 
is that the fractional average appearance in Eq.(|TP|) does not allow to use the ordinary 
Lyapunov's method because the exponential time-dependence becomes invalid. Instead, 
let us introduce the generalized exponential form: 



e^^ ^ E,{t) ^ [1 + {1 - q)tf'-'^ (24) 
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where g is a generalized Lyapunov index. First, such a type of generalization was used by 
Tsallis [8] to obtain the ordinary Gibbs-Boltzmann exponent in the limit g — > 1. In our 
case, the latter is arbitrary so that function Eq{t) acquires the power-law character, in 
particular the follow derivation rule is fulfilled: 

dE,{t)/dt = {E,{t)Y ^ El{t). (25) 

In the limits of the short and long times this function has the asymptotic behavior: 

limE,(t) ^ 1 + t, hm E,{t) ^ ((1 - q)tf'^-'' . (26) 

t— >0 t— >oo 

Below, the first of the asymptotes will be used for extraction of the Lyapunov-type mul- 
tipliers, the second one allows us to set an index q (see after Eq.(|S(]|)). 

Let us define the solutions of Eqs.(||), ( p^Qf ) in the form 

r]{t) = mE^{t), S{t) = So + nE,{t) (27) 

where 5*0 = e/3 to correspond to the point Co; the exponents /i, v and coefficients m, n 
must be determined. I Inserting Eqs. (pT]) into Eqs.(|), we obtain up to the first power of 
m, n <^ 1 



2,nEl-^'{t)E,{t) = -1. (28) 
From this, within the long-time approximation (PU]), one obtains 

n-i = 3(l-/i), u = 2. (29) 
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It is worth to note that the parameter n <C 1 to be a second term in expansion (27), whereas 
a magnitude m can be arbitrary because it stands as a first term there. PhysicaUy, it means 
that the order parameter behaves in non-Unear manner, in contrast to a linear regime of the 
structure factor. 
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Respectively, the inserting Eq.(]27|) into Eq.([T9D gives 



El-^it) [2en - aimP' EP'{t)E;\t)\ = -n. (30) 

As has been pointed out above, in the short-time hmit the function El^'^{t) can be taken 
as 1 to correspond to extraction of the Lyapunov multipher. Respectively, in the long-time 
limit there is E^^(t)E~^{t) = const = p^^ and we obtain 

3aimP' = -{l + 2e), = 1 + pi = I + [2(1 - a)]-\ (31) 

Thus, within the long time-approximation, the structure factor 

S{t) = So + {2/3){l-a)t-\ t^oo (32) 

tends to the stable magnitude So hyperbolically. The order parameter is decreased ac- 
cording to the power dependence 

r7(t) =7^0- [2(l-a)]2(i~")|m|r2(i-«), t ^ oo (33) 

with the exponent decreasing with the parameter a; the amplitude m is given by the first 
equation (PT|). 

Thus, in accordance with the phase portraits in Fig. 3a, at e < Eq the order parameter 
?7(t) decreases monotonically in the course of time, whereas the structure factor can vary 
non-monotonically, in contrast to the case a > 1/2 (cf. Fig.l). More complex behavior 
appears within the domain e > when the ordered state occurs due to the bifurcation. 
As is seen from Fig. 3b the phase plane divides into two domains corresponding to small 
and large values of the order parameter. Within the former, the system behaves as at the 
above case e < Eq, but if an initial magnitude of the order parameter is more a critical 
value, the system passes to the attractive point C. The corresponding dependencies of the 
time are depicted in Fig.6. It is characterically that these dependencies display a critical 
slowing-down near the separatrix CqSC in Fig. 3b (see curves 1, 2 in Figs. 6a, 6b). 

In order to analyze the long-time behavior in the vicinity of the point C, we can not use 
the solution like the generalized exponent (|2^). The latter is applicable when nonlinearity 
effects are sufficient to fix the above mentioned amplitudes m, n by Eqs. (p9D , (pi]). In the 
case under consideration, the linear conditions are satisfied and instead of the generalized 
exponent (0) we ought to use the Mellin transformation. The principle difference is that 
the former is defined by the single index q, whereas the later contains a set of q. Inserting 
the definitions (cf. Eqs.(p^)) 

V{t) =Vo + I m.t'^dq, (34) 
Sit) = So + I n.t'^dq (35) 
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into Eqs.(|]), (|T^) being written within the hnear approximation, we obtain the equations 
for specific amphtudes mg,nq <^ 1: 



{q/t + 2rjl)mq - 3r]ong 



0, 



[4r?o(£ 



ctiVo^ ^Pi]mq + [q/t + 2{e + r]'^)]ng 



0. 



(36) 
(37) 



This system has solutions provided the ratio c = —q/t is determined by equation 



1 ± 



\ 



1 - 



e + 2r/2 



(3J 



Thus, in a vicinity of the ordered point C the dependencies of the order parameter and 
structure factor (|33), (|35D take the forms 



(39) 
(40) 



'qit) = rjQ + m exp (— ct In t) , 
S{t) = S'o + nexp {—ct Int) , 
where the amphtudes m, n correspond to the index q = —ct. 

FinaUy, the time dependencies of the Green function G{t, 0) are determined by Eqs 



(H), (|T9D as is shown in Fig. 7. At e < sq, the monotonic decrease appears for large value of 
initial magnitude r){0) of the order parameter (see Fig. 7a). In the case e > Eq, a. maximum 
of dependence G(t, 0) disappears at condition r]{0) > rjc corresponding to the domain of 
the ordering state (see Fig. 7b). 



6 Conclusion 



Summarizing, we have derived and analized the equations for order parameter, structural 
factor and correlation function to describe the evolution of the system with multiplicative 
noise. In representation of the noise amplitude as the power- low function g{x) = x", we 
have discussed a way to apply the cumulant expansion for average {x"^"-) in the case of a 
self-similar phase space. 

We have shown that the system behavior is given by magnitude of the exponent a: at 
a > 1/2 the system is disordered; at a < 1/2 a phase transition to the ordered state 
is observed. In the first case the obtained time dependencies for the structure factor 
show the monotonic increasing to a stationary state, whereas the correlator shows the 
non- monotonic behavior. In the second case (a > 1/2) we would have to investigate the 
system behavior on the phase plane given by order parameter t] and structure factor S. It 
was shown that the phase portrait falls into two domains characterized by disordered and 
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ordered states (the former corresponds to small initial values of r]. the latter - to finite 
values of both rj and S). We have found how the system attains the steady states at long- 
time asymptotics. Within disordered domain, the structure factor decays hyperbolically 
and the order parameter dependence is described by the power-law function with the 
exponent —2(1 — a). In the ordered phase, the order parameter and structure factor 
exhibit exponential behavior with index proportional to —tint. 
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Captions 



Fig.l. Structure factor behavior at a > 1/2: 

a) time dependencies S{t) (curves 1, 2, 3 correspond to a = 0.6, e — 0.2; a — 0.6, e — 0.4; 
a = 0.9, e = 0.2); 

b) stationary point 5*0 vs temperature e for several values of exponent a. 

Fig. 2. Correlator G(t, 0) vs time at a > 1/2 for several exponent a and temperature e 
(curves 1, 2, 3 correspond to: a — 0.6, e — 0.4; a — 0.6, e — 0.2; a — 0.9, e — 0.4). 

Fig. 3. Phase portrait at a > 1/2: 

a) a = 0.3, e = 0.2; 

b) a = 0.3, e = 0.4. 

Fig. 4. Stationary states of the system at a > 1/2: 

a) order parameter r] vs temperature e for several exponent a; 

b) structure factor S vs temperature e for several exponent a. 

Fig. 5. Phase diagram eo{a). 

Fig. 6. Time dependencies corresponded to different trajectories on the phase portrait in 
Fig.3: 

a) ?7 vs Int at a = 0.3, e = 0.2 and S{0) — (curves 1, 2, 3 correspond to 7^(0) — 0.057, 
77(0) = 0.066, 7?(0) = 1.0); 

a) S ys Int a,t a — 0.3, e — 0.4 and S{Q) = (curves 1, 2, 3 correspond to 77(0) = 0.057, 
7^(0) = 0.066, rj{0) = 1.0). 

Fig. 7. Correlator G{t, 0) vs time t at: 

a) a = 0.3, £ = 0.2, S{0) = (curves 1, 2, 3 correspond to 7/(0) = 0.098, 7/(0) = 0.4, 
7/(0) = 0.97); 

b) a = 0.3, £ = 0.4, S{0) = (curves 1, 2, 3 correspond to 77(0) = 0.057, 77(0) = 0.066, 
77(0) = 1.0). 
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Fig.7 



